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Q ; ABSTRACT 

Context. B[e] stars are hot stars surrounded by circumstellar gas and dust responsible for the presence of emission lines and IR-excess 
in their spectra. How dust can be formed in this highly illuminated and diluted environment remains an open issue. 
Aims. HD 62623 is one of the very few A-type supergiants showing the B[e] phenomenon. We studied the geometry of its circum- 
stellar envelope in the mid-infrared using long-baseline interferometry, which is the only observing technique able to spatially resolve 
objects smaller than a few tens of milliarcseconds. 
Methods. We obtained nine calibrated visibility measurements between October 2006 and January 2008 using the VLTI/MIDI instru- 
(-H ' ment in SCI-PHOT mode and PRISM spectral dispersion mode with projected baselines ranging from 13 to 71 m and with various 

position angles (PA). We used geometrical models and physical modeling with a radiative transfer code to analyze these data. 
Results. The dusty circumstellar environment of HD 62623 is partially resolved by the VLTI/MIDI even with the shortest baselines. 
The environment is flattened (a/b~1.3±0.1) and can be separated into two components: a compact one whose extension grows from 
17 mas at 8pm to 30 mas at 9.6pm and stays almost constant up to 13pm, and a more extended one that is over-resolved even with the 
shortest baselines. Using the radiative transfer code MC3D, we managed to model HD 62623's circumstellar environment as a dusty 
disk with an inner radius of 3.85+0.6 AU, an inclination angle of 60±10°, and a mass of 2.10 7 M Q . 

Conclusions. It is the first time that the dusty disk inner rim of a supergiant star exhibiting the B[e] phenomenon is significantly con- 
strained. The inner gaseous envelope likely contributes up to 20% to the total N band flux and acts like a reprocessing disk. Finally, the 
hypothesis of a stellar wind deceleration by the companion's gravitational effects remains the most probable case since the bi-stability 

. . mechanism does not seem to be efficient for this star. 

. 
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(N : 

i— I ; 1 . Introduction In this context, HD 62623 (3 Puppis, HR 2996) was classi- 
cs . fled as an A2Iabe in the Yale Bright Stars Catalogue (Hoffeit & 
O ■ The B [e] phenomenon is defined by strong Balmer lines in emis- j asc hek 1982) making it one of the latest stars showing the B [e] 
^ ; sion, low excitation emission lines (especially of ionized metals phenomenon. Other classifications range from A2Ia to A3IIpe. 
K~J ' SUCh aS FeII) ' forbidden emission lines of [Fell] and [OI], and i ts distance was estimated to be 700 + 150 pc, and its radius 
^ ■ strong infrared excess due to hot circumstellar dust (Lamers et al. =55+ 10R o (Bittar et al. 2001 ). It exhibits radial velocity vari- 
H I 1998). However, these objects are not a homogeneous group of at ions with a period of 137.7 days, interpreted as direct evidence 
rt , stars in terms of stellar evolution, as such features have been de- f or an unseen companion (Lambert 1988). Plets et al. (1995) 
tected for both pre-main sequence and evolved stars. Thus, this f oun d that a 161.1-day period can also fit these variations. In 
group of stars has been divided into four sub-classes : B[e] su- both hypotheses, the system seems to be composed of a massive 
pergiants, Herbig AeB[e] stars, compact planetary nebulae B[e]- primary star with 3 1 M < M* < 39 M Q and a small mass 
type stars, and symbiotic B[e]-type stars. Moreover, since the ratio of 0.03 < M 2 /Mj < 0.15, while the projected semi-major 
spectral class of some of these stars was still very difficult to axis is 1.6 AU < a. sin i < 2.4 AU. 
determine accurately, a fifth class was added: unclassified B[e] 

stars. Recently, in an effort to investigate this mysterious sub- . T ^ Ts ,,~,„x ,. 

group, Miroshnichenko (2007) proposed to classify most of their u Using .color indices (H - K) and (K - L), Allen ( 973) found 

members as a new homogeneous group of non-supergiant and * at H ? 6 } 62 l s lnfrared ex = ess f t0 ° lar § e l ° onl y orl g lnate 
non-pre-main sequence stars. He proposed to name this group from the free - free and free-bound emissions of a gaseous en- 
after its precursor, i.e. FS CMa's group. Finally, the binarity de- ^ elo P e and 18 P"*aMy due to the presence of optically thick 
tected for 30% of these stars could play a key role for the origin dust ln the "rcumstellar environment of this star. Observations 



of the B[e] phenomenon for these FS CMa stars. ' with the Low-Resolution Spectrograph (LRS) of the InfraRed 

Astronomical Satellite (IRAS) clearly exhibit a 9.6 /im silicate 
feature in the object spectra, confirming the presence of dust in 



Send offprint requests to: meilland@mpifr-bonn.mpg.de the envelope. Since the star also shows evidence of Ha and [OI] 

* Based on observations made with ESO Telescopes at Paranal emission lines (Swings & Allen 1976), it was classified as a B[e] 
Observatory under programs 078.D-O51 1 and 080.D.0181. supergiant, i.e. sgB[e] (Miroshnichenko 2007). 
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Fig. 1. Left : HD 62623 SED from various sources of literature: red diamonds are flux measurements listed in Bittar et al. (2001); 
green lines, IRAS LRS spectra; dark blue lines, IUE spectra; light blue lines with errors bars, our VLTI/MIDI measurements; and 
the dark line is the best reddened Kurucz model, whose parameters are presented in Table 1 . Right : Zoom of the N band of the SED 
(8-13^/m) to exhibit differences between the various measurements obtained in this spectral band and the VLTI/MIDI spectra. 



The formation mechanisms responsible for sgB[e] stars' cir- 
cumstellar environment are still an open issue. Spectroscopic 
and polarimetric observations suggest that their circumstellar 
envelopes are non-spherical (Magalhaes 1992). Zickgraf et al. 
(1985) proposed a general scheme for these stars consisting of 
a hot and fast line-driven polar wind responsible for the pres- 
ence of forbidden lines and a slowly expanding dense equatorial 
region where permitted-emission lines and dust can form. The 
deviation from the spherical geometry for this model can be re- 
produced by a nearly critically rotating star. 

Using aperture masking interferometry with the Keck I tele- 
scope, Bittar et al. (2001) partially resolved HD 62623's cir- 
cumstellar envelope between 1.65 and 3.08 fim and showed that 
the envelope extension was increasing from 6-7 mas at 1.65 
fim to 25-30 mas at 3.08 fim. Stee et al. (2004) have success- 
fully modeled both these observations and the Spectral Energy 
Distribution (SED) using a gas + dust model with an inner dust 
radius of 15-20R* and a temperature at this boundary of 1500K. 
In their best model, the gas emission originates from the inner 
part of the envelope and not from a polar wind. Nevertheless, 
the issue of HD 62623's envelope flattening remains unsolved, 
even if Yudin & Evans (1998) measured a polarization in the vis- 
ible of 1.5% with a polarization angle of 95+5°, which seems to 
favor a non-spherical object. 

By putting strong constraints on the envelope geometry, long 
baseline infrared interferometry can be a key technique to study 
the physics of active hot stars (Stee et al. 2005). Using both 
VLTI/MIDI and VLTI/AMBER data, Domiciano de Souza et 
al. (2007) resolved the circumstellar environment of the sgB[e] 
CPD -57° 2874 and modeled it as a flattened compact envelope 
with an extension increasing with wavelength. 

In this paper, we present new spectrally resolved VLTI/MIDI 
interferometric observations of HD 62623 that allow us to char- 
acterize its circumstellar envelope geometry and determine its 
physical parameters. 

The paper is organized as follows. In Section 2 we start our 
modeling of HD 62623 by reconstructing its SED in order to 
constrain its stellar parameters and distance. The VLTI/MIDI ob- 
servations and the data reduction process are briefly introduced 
in Section 3, and a first analysis of the calibrated visibilities us- 
ing simple geometrical models is presented in Sect. 4. The result 
of these study is used as a starting point for a three-dimensional 



Monte-Carlo radiative transfer modeling presented in Sect. 5. 
Finally, the physical properties of HD 62623 are discussed in 
Sect. 6, while Sect. 7 draws the main conclusions of this study. 

2. Preliminary study 

2. 1 . The Spectral Energy Distribution 

In order to accurately determine the stellar parameters of 
HD 62623, we started our study by reconstructing the ob- 
ject SED using photometric and low resolution spectral mea- 
surements ranging from ultra-violet to far-infrared. We used 
data from various sources in the literature already described in 
Bittar et al. (2001). These data include: IRAS fluxes and Low 
Resolution Spectra (LRS), IUE spectra, and other photometric 
measurements in the UV, the visible, and the near-IR. We also in- 
cluded our N band VLTI/MIDI spectra in this study. The result- 
ing reconstructed SED is plotted in Fig.Q] The strong IR-excess 
due to circumstellar dust is evidenced by the inflection point in 
the SED around 1 .5 fim and by the 9.6 //m silicate emission band 
visible on IRAS LRS and VLTI/MIDI spectra between 8 and 13 
fim. The VLTI/MIDI flux is of the order of the other mid-IR mea- 
surements. This clearly shows that most of the N-band emission 
originates from the 250 mas interferometric field of view of this 
instrument. We note that a second silicate emission band is also 
visible in IRAS LRS spectra around 20 fim. 

2.2. Circumstellar or interstellar extinction 

We have tried to fit the UV and visible parts of the reconstructed 
SED using typical Kurucz models (Kurucz 1979) for supergiant 
stars; i.e, with 1.5 < log g < 2.5 and 8000K < T eff < 10000K. 
None of these models were able to simultaneously fit all the SED 
from UV wavelengths up to 1 fim, and the observed SED seems 
to be slightly reddened. To explain this reddening, two hypothe- 
ses can be formulated: (1) if the circumstellar dust becomes opti- 
cally thick at short wavelengths, it can absorb a part of the stellar 
UV emission, assuming this dust is located in the line of sight, 
or (2) the reddening can be due to interstellar matter. 

It is hard to discriminate between these two hypotheses 
without fully modeling HD 62623's circumstellar environment. 
However, it is still possible to estimate the reddening without 
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Villus 




8250 ± 250 K 




2.0 ± 0.5 


R* 


65 ±5Ro 


d 


650 ± 100 pc 


E(B-V) 


0.17 ± 0.03 (1) 



Table 1. Stellar parameters, distance, and extinction determined 
from the fit of HD 62623 SED. (1) value from Plets et al. (1995) 



solving the issue concerning its origin. Using Meynet & Hauck 
(1985) intrinsic Geneva colors laws for A-F supergiants, Plets et 
al. (1995) estimated a reddening E(B-V) of 0.17+0.03. 



2.3. Stellar parameters 

Assuming that the extinction originates from a standard interstel- 
lar medium (Ry=3.1) and that it follows Cardelli et al.'s (1989) 
law from UV to far-IR wavelengths, the SED can be success- 
fully fitted with reddened Kurucz models with T e ff ranging from 
8000Kto 8500K and log g from 1.5 to 2.5. Our best model, using 
T e yy=8250K and log g=2.0, is overplotted in Fig. [TJ The good 
agreement between the observed SED and the modeled one is 
obtained for an c//R* ratio of 10+1, where the distance d is ex- 
pressed in parsecs and the stellar radius R* in R Q . Since a typical 
value of R* for a A supergiant is 65 ±5 (Allen 2000), we can es- 
timate a distance for the star of 650± 1 00 pc, which is compatible 
with the previous estimation by Bittar et al. (2001). The resulting 
parameters used in this paper to model HD 62623 's central star 
are presented in Table 1 . 



3. VLTI/MIDI Observations and data reduction 

VLTI/MIDI observations of HD 62623 were carried out at 
Paranal Observatory between October 2006 and January 2008 
(Guaranteed Time Observing runs 078.D-051 1 and 080.D.0181) 
with the 1 .8m Auxiliary Telescopes (ATs). We obtained nine vis- 
ibility measurements with projected baselines ranging from 13.4 
to 71.4 meters and with various orientations on the sky plane. 
Three stars were used as calibrators during this observing cam- 
paign: Procyon, Sirius, and Alphard. The log of these observa- 
tions is presented in Table 2. 

All observations were made using the SCI-PHOT mode, 
which enables a better visibility calibration since the photom- 
etry and the interferometric fringes are recorded simultaneously. 
Thanks to the PRISM low spectral dispersion mode, we also ob- 
tained spectrally resolved visibility with R=30 in the N band 
(7.5-13. 5/zm). Two different reduction packages were used to 
reduce these data : MIA developed at the Max-Planck Institut 
fur Astronomie and EWS developed at the Leiden Observatory 
(MIA + EWS, ver. 1.5.1). Since the two methods provide the 
same results within the error bars, we decided to use the MIA 
reduced visibilities in the following work. 

HD 62623 's (u,v) plane coverage is plotted in Fig [2] The 
large spread of baseline length (between 13.4 and 71.4) and PA. 
(between -4.1° and 85.3°) will allow us to put strong constraints 
on the envelope geometry in the N band (i.e. extension, flatten- 
ing, radial intensity profile). Moreover, two triplets of baselines 
are almost aligned: B 2 , B 3 , and B 6 with PA.~40°, and B 4 , B 5 , 
and B9 with PA.~85°. This will simplify our modeling by allow- 
ing us to directly determine the object projected intensity profile 
at these two orientations. 
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Fig. 2. UV plane coverage obtained for HD62623. The dot- 
ted line represents the polarization angle measured by Yudin 
& Evans (1998). The circle colors represent the two triplet of 
aligned baselines discussed in Sect. 3 : red for PA~40° and green 
for PA~85°. 
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Fig. 3. HD 62623 VLTI/MIDI calibrated visibilities for the nine 
baselines from Fig. 2 plotted as a function of wavelength. The 
dotted line corresponds to the central wavelength of the silicate 
band, i.e, 9.6 /im. 
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Name 


wds. uate 


Obs. Time 
(U 1) 


Stations 


Projected baseline 
Length(m) P.A.(°) 


Calibrators 


Bi 


2006/ 10/1 / 


06:42 


D0-H0 


49.5 


27.2 


7 Eri, Procyon, /? Gru 


B2 


2006/10/1 / 


08:10 


D0-H0 


56.8 


/ion 
48.0 


7 Eri, Procyon, p Gru 


B 3 


2006/12/16 


04:00 


G0-H0 


27.8 


45.8 


Procyon, 4 Eri 




zuuo/ iz/ ly 


Uo . J 




56.6 


84.8 


Alphard, Procyon 


Bs 


2006/12/21 


08:19 


G0-H0 


29.0 


84.0 


7 Eri, Procyon, Alphard 


B 6 


2007/10/04 


08:22 


E0-G0 


13.3 


39.4 


Sirius 


B 7 


2007/12/10 


05:22 


G1-H0 


71.4 


-4.1 


Sirius 


B« 


2007/12/14 


08:13 


G1-H0 


70.9 


15.4 


Sirius 


B 9 


2008/01/12 


07:08 


E0-G0 


14.0 


85.3 


Sirius, Alphard 



Table 2. VLTI/MIDI Observations log of HD 62623 and its calibrators. 



4. Results 

4. 1 . Visibilty and Gaussian equivalent FWHM 

The calibrated visibilities for the nine projected baselines are 
plotted as a function of wavelength in Fig. [3] All these spec- 
trally resolved visibilities, except the ones from the two longest 
baselines (i.e., B7 and Bg), exhibit a small drop with a minimum 
around 9.6 +0.1 /mi. Such an effect, often found in N-band in- 
terferometric observations of dusty circumstellar environments, 
is due to an opacity increase between the continuum and a sil- 
icate band. This greater opacity in the silicate band causes the 
circumstellar envelope to appear more extended than in the con- 
tinuum. The fact that such an effect is not visible in the B7 and 
Bs data probably means that both the silicate and the continuum 
emission are over-resolved with these baselines. We note that, 
for this A[e] supergiant star, the silicate band is clearly in emis- 
sion (see Fig.[TJ. 

The object is partially resolved for all baselines, even for 
the shortest ones (i.e., V<0.9 for B 6 =13.3m and B 9 =14.0 m). 
This implies that the circumstellar envelope is quite extended. 
To roughly estimate its size, we calculate the Gaussian equiva- 
lent FWHM for each baseline and at three wavelengths (i.e., 8, 
10, and 12 /zm). The results, presented in Table 3, clearly show 
that the envelope extension is growing by about 45% between 8 
and 10 /mi, whereas it remains nearly constant between 10 and 
12 /mi. As we already mentioned, the envelope size variation 
between 8 and 10 /mi is mainly due to the change of opacity 
between the continuum and the silicate band, but another effect 
has to be taken into account to explain why the size does not 
drop back to its 8 /mi value from 10 to 12 /mi. Remembering 
that larger wavelengths can probe cooler region it implies, for 
an untruncated gaseous or dusty envelope heated by a central 
source, i.e., a single star, that the size of the emission is grow- 
ing with the wavelength as the cooler regions are located further 
from the heating source. If the size in the continuum were not 
increasing with the wavelength this would have implied that the 
disk was somehow truncated, as it was discovered for the clas- 
sical Be star a Arae (Meilland et al. 2007), or that the disk was 
mainly isothermal (Jones et al. 2004). 

If the intensity distribution of the object was close to an el- 
liptical Gaussian, the estimated FWMH should only depend on 
the baseline orientation and not on its length. As we already 
mentioned in the previous section, two groups of three base- 
lines are almost aligned. Consequently, the Gaussian equivalent 
FWHM should be nearly constant within each one. As is shown 
in Table 3, this is not the case for HD 62623 and for a given 
wavelength and P. A. The estimated size shrinks with the baseline 
length and finally, the Gaussian equivalent size measured for the 
shortest baseline is about twice the one measured for the longest 
baseline. This result is also valid if we model the object as a sin- 





Gaussian Disk FWHM (mas) 


P.A.(") 




8//m 


lO/zm 


12/mi 




B 4 


17.6+ 0.3 


24. 1± 0.3 


24.4± 1.2 




B 5 


21.5+ 0.7 


30.0± 1.5 


32.3± 1.0 


P.A.~85° 


B 9 


31.7+ 4.0 


40.7± 4.9 


43.4± 6.3 




B 2 


16.6± 2.6 


26.4± 1.0 


26.1± 1.5 




B 3 


22.7± 0.7 


32.8± 1.7 


31.6+ 3.5 


PA.~40° 


B 6 


43.8± 3.8 


58.7+ 4.8 


59.2± 5.9 




Bi 


21.2+ 2.6 


32.4± 1.1 


32.4± 1.6 


P.A.=27.2° 


B 7 


16.8+ 0.7 


23. 8± 0.7 


27.7± 1.0 


p.a.=-4.t 


B 8 


15.4± 0.6 


22.4± 0.6 


24.9± 0.7 


P. A. = 15.4° 



Table 3. Gaussian disk equivalent FWHM for all baselines and 
at three different wavelengths. The baselines are classified by 
P.A. in order to exhibit the inconsistency of a single Gaussian 
modeling, which would have provided similar sizes for baselines 
with the same P.A. 
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Fig. 4. HD 62623 VLTI/MIDI calibrated visibilities plotted as a 
function of the spatial frequency B/A. The colors correspond to 
the orientation of the projected baselines: green for 85°, blue for 
40°, and red for 27.2°. 17 mas (dashed line) and 60 mas (dotted 
line) Gaussian disks are overplotted to illustrate the inconsis- 
tency of modeling HD 62623's circumstellar environment as a 
single Gaussian disk or ellipse. 
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Fig. 5. Best-fit models parameters and^- 2 for a fully resolved component + elliptical Gaussian geometry (Gl). From left to right: 
Normalized flux of the most compact component Fi, FWHM of the major and minor axis of the compact component, and corre- 
sponding reduced^ 2 . The vertical dotted line indicates the position of the center of the silicate band, i.e 9.6 yum. 



gle uniform ellipse or elliptical ring. This clearly shows that the 
object geometry is more complicated and should be modeled us- 
ing at least two components: a compact one which is resolved 
only by the longest baselines and a more extended one which 
may even be fully resolved by the smallest baselines. 

This assumption becomes obvious when plotting the visibil- 
ity for all baselines as a function of the spatial frequency \B\/A 
(Fig|4]i. This plot also exhibits all the previously described char- 
acteristics of the object in a very synthetic way. The wavelength 
dependence is evidenced in each of the seven curves representing 
the visibility for each baseline. These curves are inverted com- 
pared to Fig|3](i.e., the wavelength increases from right to left), 
but both the silicate band and the continuum effects on the visi- 
bility are shown. Clues of the flattening of the envelope, at least 
in the silicate band, are also evidenced since the visibility de- 
pends not only on the baseline length but also on its orientation. 
The object is larger in the B2, B3, and B9 orientation (P.A.~40°) 
than in the B4, B5, and B9 one (P.A.~85°), and is even larger 
along Bi (PA. =27.2°). These visibility variations as a function 
of the P.A. are consistent with an orientation of the major-axis 
of 5° derived from Yudin & Evans (1998) polarization measure- 
ments. Nevertheless, the flattening cannot be accurately deter- 
mined without knowing the radial intensity profile of the object. 
We note that B7 and Bg are almost aligned with the putative en- 
velope major-axis. 

4.2. A two-component model 

It is obvious from Table 3 and Figg]that HD 62623 MIDI visi- 
bilities cannot be fitted by a single component model. Thus, we 
try to fit them with three different kinds of two-component geo- 
metrical models. For all of these models, the most compact com- 
ponent is an elliptical Gaussian defined by three parameters: the 
major-axis FWHM (ai), the minor-axis FWHM (bi), and the ori- 
entation of the major-axis into the sky plane (9\). Since the po- 
larization angle measurement seems to provide the orientation 
of the minor-axis with a good accuracy, we can set the value of 
61 to 5°. We will see in Sect. 5 that this component probably cor- 
responds to the inner rim of the circumstellar dusty disk and that 
it can be modeled as an elliptical ring. However, the envelope is 
not fully resolved and this component can be successfully fitted 
by several other flattened models including Gaussian or uniform 
ellipse. 

The last parameter concerning this compact component is the 
flux ratio between this Gaussian (Fj) and the total normalized 
flux at each wavelength. The flux of the extended component 



is consequently (1-Fi). Since we want to determine if it can be 
fully resolved by the smallest baselines, we model it following 
three different ways of increasing complexity to test if we can put 
some constraints on its geometry. To check the consistency of 
our two-component modeling, we also add a "reference" model 
assuming only one elliptical Gaussian. The parameters of the 
four different geometries tested (GO, Gl, G2, and G3) are the 
following : 

- GO : This is the "reference" model. The geometry is defined 
by a single elliptical Gaussian. This model has only 2 free 
parameters : ai, and b\. 

- Gl : Here, we consider that the extended component is fully 
resolved, and we model it as a homogeneous constant con- 
tribution. This first model has 3 free parameters : ai, bi, and 
Fi. 

- G2 : We try to model the extended component as a Gaussian 
with the same flattening as the compact one so that it is de- 
fined only by its major axis FWHM a2. This model has 4 free 
parameters : ai, &2> bi, andFi. 

- G3 : In this last model, the minor-axis FWHM (b2) of the ex- 
tended component is also a free parameter. Thus, this model 
has 5 free parameters : ai, aa, bi, b2 and Fi. 

Since the envelope extension varies within the N band, we 
assume that the free-parameters are all wavelength dependent. 
For each wavelength, we try to minimize the model reduced x 1 
using a Levenberg-Marquardt 2D algorithm with 10 3 random 
starting positions. We finally calculate the global^ 2 for all wave- 
lengths to compare the best-fitted models of the three geometries 
(and the reference GO). The value of the global^ 2 are 9.14 for 
GO, 2.62 for Gl, 2.62 for G2, and 3.26 for G3. Thus, this mod- 
eling clearly shows the lack of data to constrain the geometry of 
the large component of our model but the consistency of a two- 
component modeling approach. More proof of the difficulty to 
constrain the large component geometry is the fact that the wave- 
length dependence of a2 and b2 for the best G2 and G3 models is 
erratic, with values randomly distributed between 100 mas and 
1000 mas. Finally, it is obvious that we need more measurements 
with smaller spatial frequencies to constrain the geometry of this 
extended component. Nevertheless, we can still determine lower 
and upper limits to its size considering that it is fully resolved 
at 8 //m with the 13 m baseline (~150 mas for a uniform disk 
model) and smaller than the VLTI/MIDI interferometric field of 
view (~250 mas at 10/mi). 

On the other side, the relative flux of the two components 
and the minor and major axis of the compact component are 
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well constrained and do not depend on the geometry used to 
model the extended component in the 8-12.5 pm range. We have 
shown that we cannot constrain the geometry of the extended 
component; consequently, in the following, we will discuss only 
results concerning the best fitted Gl model. The wavelength de- 
pendence of Fi, ai, ai/bi, and the reduced^ 2 for the best-fitted 
Gl model are plotted in Fig. [5] 

This figure clearly exhibits the effect of the silicate band cen- 
tered around 9.6 pm on the object geometry. The compact com- 
ponent extension and flattening strongly varies between 8 and 10 
pm. At 8 pm this structure does not show any flattening and ai 
and bi are about 15-16 mas. Their respective size increases to 
about 22 and 28 mas at 10 pm and then remains almost constant 
up to 13 pm. On the other hand, the relative flux of the compact 
component Fi/F for is almost constant between 8 and 10 pm and 
of the order of 0.75 ±0.02. Beyond 10 pm it starts rising up to 
0.82 ±0.03 around 13 pm. 

As we already mentioned in the previous sections, the enve- 
lope extension variation in the N band can easily be explained 
by a silicate band and continuum effect. However, the apparent 
variation of its flattening remains hard to explain in the scheme 
of a classical dusty disk. If the disk is Keplerian and stratified, 
we would expect it to be flared. Remembering that larger wave- 
lengths are probing cooler regions, and thus further from the cen- 
tral star, we would expect that the apparent flattening of a flared 
disk decreases with wavelength. Another possibility is that the 
disk might be created by an outflowing radiative wind. In this 
case, if the terminal velocity of the wind is high enough to avoid 
any stratification of the disk, we would expect that the flattening 
does not depend on the wavelength. Additional effects have to 
be considered to allow the disk to appear more flattened at larger 
distances from the central star. One solution might be the pres- 
ence of a puffed-up inner rim, as already proposed for Herbig 
stars by Natta et al. (2001). However more measurements are 
needed to solve these issues. 

If we consider that the 1 3 pm envelope is a geometrically thin 
equatorial disk, the flattening would only be due to its projection 
onto the sky plane, and we can estimate HD 62623's inclina- 
tion angle to be 38±8" (using the I3fim value for the flattening). 
Unfortunately, the disk may not be geometrically thin and may 
also be flared, as we have already discussed. Consequently, we 
can only estimate a lower limit on this parameter of 30°. 



4.3. Comparison with KECK measurements 

Using Keck I aperture masking interferometry, Bittar et al. 
(2001) partially resolved HD 62623's circumstellar environment 
in the near-infrared. They estimated the circumstellar environ- 
ment extension using Gaussian disk models and found that its 
size was 6.8 ±1 mas at 1.65 fim, 24.8 ±1 mas at 2.26, and 27.3 
±3 mas at 3.08 fim. Even if their estimation at 3.08 is almost 
two times larger than our compact component at 8 fim, one must 
keep in mind that they only obtained measurements for spatial 
frequencies less than 0.02 cycles/mas, and that we should com- 
pare it to a "single" Gaussian estimation for our small and inter- 
mediate baselines (i.e., B3, B5, Be, and B9 from Table 3). Using 
these measurements, we obtain an average size of 29.9 mas at 
8 fim. This value is clearly compatible with the one determined 
by Bittar et al. (2001) at 2.26 and 3.08 fim. Moreover, the esti- 
mated size does not strongly vary between 2.26 and 8 fim (i.e., 
only 1.21 times larger), which might indicate that the measured 
extension in the K band is close to the diameter of the inner rim 
of the dusty disk. 




Fig. 6. Typical radial and vertical density distribution for an 
equatorial wind (left) and a Keplerian viscous disk (right). 

The existence of the two components in our HD 62623 mod- 
eling is probably due to the presence of a bright inner rim in 
the dusty disk which probably lies within the compact compo- 
nent of our model and is "blurred" as we do not fully resolved it. 
Thus, the inner radius of the dusty disc should be smaller than 
the Gaussian equivalent size measurement of the compact com- 
ponent at 8 fim. 

We obtained 15 mas ~ 16 R+ at 8 fim for the compact com- 
ponent Gaussian equivalent FWHM, and the extension at 2.26 
fim is 1.21 times smaller than the one at 8 fim. Moreover, when 
the envelope is not fully resolved, estimating its size with a ring 
diameter corresponds to a value 1.24 times smaller than if we 
were using a Gaussian FWHM. Thus, we can roughly estimate 
the dusty disk inner rim radius to be of the order of R,„ ~10R* 
or 3 AU. 

5. MC3D simulations of HD 62623 

5.1. The Monte-Carlo code 

MC3D (Wolf et al. 2003) is a three-dimensional continuum ra- 
diative transfer code. It is based on the Monte-Carlo method and 
solves the radiative transfer equations self-consistently. Starting 
with a spatial density distribution of scatterers and absorbers and 
some primary radiative source(s), it can be used to determine 
the temperature in this medium by an iterative process of heat- 
ing, absorption, and reemission by dust grains. When the self- 
consistent temperature distribution is reached, the software can 
compute the resulting observables: SED, wavelength-dependent 
images, and polarization maps. 

The software allows various geometries for the density dis- 
tribution : one-dimensional (i.e., dust shells) , two-dimensional 
(i.e, axi-symmetric disk or torus), or fully three-dimensional in 
spherical, cylindrical, or Cartesian coordinates (for example, to 
model inhomogeneities in a circumstellar disk). 

5.2. Dust disk model 

We decided to model HD 62623 with a two-dimensional dust 
density distribution and a radial and vertical structure, as de- 
scribed in Shakura & Sunyaev (1973) and Wood et al. (2002) : 



p(r, z) = po 



exp 



2h(r) 2 J 



(1) 



where r is the radial distance in the midplane of the disk, a 
the density parameter in the midplane, and h(r) is given by : 



h(r) = h 



R* 



(2) 
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Fig. 7. From left to right : 8.2, 9.9, and 13.1 //m MC3D intensity maps for the best-fitted equatorial wind model (top) and Keplerian 
disk (bottom) of HD62623's circumstellar environment. 



where /3 is the vertical density parameter, and ho is the theo- 
retical disk scale height at R* . In order to be more understand- 
able, the scale height at 100 AU (hioo) is used as an MC3D input 
parameter instead of ho. The three geometrical parameters a, /?, 
and ho are related to the physical nature of the dusty disk. Porter 
(2003) shows that both the viscous disk and the equatorial wind 
model can theoretically allow the formation of dust around su- 
pergiant B[e] stars. A Keplerian viscous disk will imply a=3.5, 
/3-1.5, and ho=c s / VGM*, where c s is the sound speed. On the 
other hand, a -2 and (3-1 for a wind model with a constant ex- 
pansion velocity. In this case, ho roughly defines the opening an- 
gle of the dusty disk. A typical density distribution for an equa- 
torial wind and a Keplerian viscous disk are plotted in Fig [6] 

Details of the spatial distribution of chemical compositions, 
shapes, sizes, and other parameters of real dust grains in this 
object are poorly known. Large uncertainties of the real proper- 
ties of dust enormously increase the free parameter space of the 
models to be explored. Since our purpose is to reconstruct global 
properties of the disk, we decided to fix the dust grain properties. 
We have only treated the transfer of dust radiation in this work. 
The gas component, present in the model only implicitly, is de- 
scribed by a dust-to-gas ratio of 1 % within the dense disk. 

We assumed the standard interstellar size distribution from 
Mathis etal. (1977): 



3n(a) _ 3 5 
— — a a 
oa 



(3) 



where a is the dust grain radius extending from 0.005 to 
0.25/mi (Draine & Lee 1984). We considered the dust grains to 
be homogeneous spheres. 



5.3. Model fitting 

We calculated several hundreds of models to explore the free- 
parameters space and we tried to constrain their values and un- 
certainties. We focused on both the viscous disk and the equato- 
rial wind models defined in the previous section to test if we can 
clearly discriminate between them. Finally, we also tried many 
intermediate models which imply that the disk is slowly expand- 
ing and not perfectly stratified. 

We considered the following outputs from our simulations: 

1. The total SED computed from 1 to 20 fim that can be di- 
rectly compared to our reconstructed SED, including our 
MIDI data presented in Sect. 2. 

2. Twenty 301x301 pixels intensity maps in the MIDI wave- 
length range, i.e, 8-13.5 fim, with a pixel size of 0.7 mas. 
Each image allows us to fit the spectrally dispersed visibil- 
ities at one wavelength. The transformation of the images 
into a visibility signal is straightforward: the images are ro- 
tated by the PA angle of the object onto the sky and then 
collapsed as ID flux distribution in the direction of each of 
the nine baselines. These vectors are Fourier transformed and 
normalized. 



For each model we compute the SED and visibilities x 2 sep- 
arately. The use of a single global x 1 would have implied solv- 
ing the problem of the relative weight of the SED and visibility 
in the computation. Finally, we decided to concentrate our ef- 
forts on the visibility fit since it contains more information on 
HD 62623's circumstellar environment geometry than the SED. 
Nevertheless, for each model we verify that the fit of the SED 
is accurate enough, both in term of IR-excess and silicate band 
emission intensity. 
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Fig. 8. Fit of the nine visibility measurements (top) and SED (bottom) of the best equatorial wind model (left) and Keplerian viscous 
disk model (right). The dotted lines correspond to the models and the errors bars to the data. 



5.4. The best models 

The parameters of the best equatorial wind and Keplerian vis- 
cous disk models are presented in Table |4] and the correspond- 
ing 8.2, 9.9, and 13.1/im images for both geometries are plotted 
in Fig [7] The fit of the visibilities is plotted in Fig [8] These fig- 
ures show that both scenarios are able to reproduce the general 
trends of the seven shorter baseline visibility variations, but the 
ones measured for B7 and Bg are significantly underestimated for 
wavelengths shorter than 1 1 /im. Such a residual visibility might 



Parameter 


Best Wind Model 


Best Keplerian disk model 




4±0.5AU 


3.7+0.5AU 


Ron? 


150AU 


150AU 


Menv 


2.10- 7 Mo 


2.1(T 7 M 


Tdust 


1250K 


1250K 


a 


2 


3.5 


p 


1 


1.5 


h-lOOAU 


24AU 


28AU 


Inclination angle 


60±10° 


60±10° 


Major-axis P.A. 


10±10° 


15±10° 



Table 4. Parameters of the best fitted MC3D equatorial wind 
and Keplerian disk model. For both models T e ff=8250K and 
L + =17OOOL (i.e., corresponding to R*=65R Q ) and d=650pc. 



suggest the presence of structures not fully resolved with the 
70 m baselines, which are not taken into account in our MC3D 
model. Some hypotheses on the nature of these structures will 
be discussed in Sect. 6.1. 

The reduced x 1 on the interferometric data of the best equa- 
torial wind and Keplerian disk are 14.9 and 11.5, respectively. 
These high values are mainly due to the poor fit of the 70m base- 
line. Thus, skipping B 7 and Bg measurements, the reduced x 1 
reach values of 2.8 and 2.4, respectively. The main qualitative 
difference between the two scenarios is the silicate band visibil- 
ity drop scale. This favors the hypothesis of the stratification of 
the disk. 

However, the differences are small, and a large set of inter- 
mediate values of a and /3 are also possible so that the disk strat- 
ification cannot be inferred accurately. Actually, as can be seen 
in Fig|7] the geometrical differences between these scenarios are 
too small to be properly constrained by our nine VLTI/MIDI 
measurements with uncertainties of the order of 5%. 

Nevertheless, despite our inability to deduce accurate val- 
ues of the radial and vertical disk structure from the visibility 
fits, we clearly constrain the inner disk radius, R,„=3.85±0.7 (or 
R in =4±0.5 AU for the equatorial wind and R,„=3.7 +0.5 AU for 
the stratified disk). This represents the first measurement of the 
extension of the inner rim of a B[e] supergiant dusty disk. The 
mass of the envelope does not strongly depend on the disk strati- 
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Fig. 9. Left : Visibility plotted as a function of the spatial frequency for several random clump distributions (colored solid lines) with 
N c iump= 10000 and FWHM e /„ mp =3 mas, using the 8 fim intensity map from our best MC3D wind model as an input. The visibility 
corresponding to this map is also plotted as a dashed line for comparison. The vertical dotted lines represent the position of the 
spatial frequencies for a 70 m baseline at 8 and 13 fim. Right : Example of a clumped map with a random distribution for the same 
parameters. 



fication and is also well constrained, i.e., 2 ±1 10~ 7 M , so is the 
inclination angle, i.e., 60 ±10°. Finally, the position of the major- 
axis that is roughly perpendicular to Yudin & Evans (1998) po- 
larization angle measurement 95°. 

The SED fit for each scenario is presented in Fig [8] For 
both scenarios, the fit of the SED is not good enough, especially 
the 9.6/im silicate peak, which is too strong for both models, 
whereas the IR-excess is overestimated for the equatorial wind 
and underestimated for the Keplerian disk scenario. However, 
the general agreement is better for the Keplerian disk, and we 
note that both models managed to reproduce the 13-20 fim SED 
well. 

In fact, an accurate fit of the silicate feature is very often dif- 
ficult to obtain. The strength of the feature depends on the size 
and composition of the silicate grains, but also on the level of the 
underlying continuum, thus, on the quality of the model density 
structure. The near-IR region of the SED is also strongly influ- 
enced by the detailed structure of the inner rim of the disk which, 
in this study, is assumed to be a wall. Recently, many studies 
were devoted to this point in the field of young stellar objects, 
involving a puffed-up inner rim and dust settling (Dullemond et 
al. 2001, 2004; Tannirkulam et al. 2007). 

For the Keplerian disk, the near-IR flux is underestimated 
because the inner rim is probably not large enough. On the con- 
trary, the equatorial wind model overestimates the near-IR as a 
direct consequence of the smaller density coefficient a, making 
the inner rim almost as opaque as a true wall. Moreover, the 
flaring of the density structure is also larger; thus, more of the 
surface of the dust disk is directly exposed to the stellar radi- 
ation and the vertical extent of the structure is probably over- 
estimated. 



6. Discussions 

6.1. Nature of the residual visibility 

As already mentioned in Sect. 4.1, the fact that the two largest 
baselines, i.e. B7 and Bg, do not exhibit the same wavelength de- 
pendence as the shorter ones, in particular the lack of a clear drop 



of the visibility in the 9.6 fim silicate band, might be a clue that 
both the continuum and the silicate emissions are over-resolved 
with these baselines. Since the visibility does not reach zero with 
these baselines but 0.17+0.05 at 8 fim and 0.12 +0.03 at 13 fim, 
this might advocate for the presence of a unresolved structure 
which contributes to the total flux at least to 17% at 8 fim and 
12% at 13 fim, assuming it is fully unresolved. 

This assumption is confirmed by our MC3D models that 
clearly show the impossibility to reproduce the B7 and Bg mea- 
surements simultaneously to the other ones. Actually, with these 
70 m baselines, the dusty disk inner rim (i.e., R; n = 3.85 AU = 
12. 8R*) is almost fully resolved. Thus, the visibility should be 
close to zero. Taking a smaller R; n would lead to a better fit of B7 
and Bg visibilities, but to an overestimation for the other ones. 
Moreover, the fit of the largest baselines would not be totally sat- 
isfying since the silicate band effect would be clearly visible in 
the modeling. 

Thus, in the following, we investigate two different hypothe- 
ses to explain the residual visibility measured with the 70 m 
baselines. 



6.1 .1 . Clumping of the dust disk 

Effects of clumping on the visibility are not trivial and depend 
mainly on the number of clumps and their extension. If they are 
unresolved, i.e. D c ; 1(mp <K B/A, they will significantly increase 
the visibility level. However, if their number is small, they will 
act like multiple objects and produce oscillations structures in 
the visibility function as a function of the spatial frequency. 
Consequently, since a clear increase of the visibility level but 
no oscillations are evident in our data, we can assume that the 
clumps are unresolved and numerous. 

To simulate the effect of clumping on our visibility mea- 
surements, we use a simple and ad-hoc method to produce 
a clumped, dusty disk emission map. We compute a random 
clumps map in which each clump is a normalized Gaussian. We 
then multiply this map by an MC3D emission map. The result- 
ing map consists of a clumped version of our previous MC3D 
emission map. Fig [9] (left) presents the visibility plotted as a 
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Fig. 10. Left : Visibility plotted as a function of the spatial frequency for the SIMECA + MC3D model with parameters from Table[5] 
Right : corresponding 8/mi SIMECA intensity map. 



function of the spatial frequency for several random clump dis- 
tributions using a 13.1 pm intensity map from our best MC3D 
Keplerian disk model as an input and with N c i ump = 10000 and 
FWHM c i ump = 3 mas. An illustration of a randomly clumped 
map is also plotted in Fig [9] (right). 

This figure exhibits the random clump distribution effect 
on the visibility. For the selected clump parameters, it clearly 
shows that when the global disk structure is fully resolved (i.e. 
B/A >0.05), the visibility can take any value between and 0.2, 
depending on the clump distribution. Such random effect due to 
dumpiness have already been studied using fully radiative trans- 
fer clumped-disk models by Honig et al. (2006) for AGN tori. 

Finally, it seems that random clump distributions with the 
previous parameters can produce visibilities compatible with our 
measurements. We note that this result is strongly dependent on 
the baseline P. A., and two baselines separated by a few degrees 
can produce totally different visibility profiles. Consequently, 
knowing that our two largest baselines are separated by nearly 
20°, the fact that they exhibit nearly the same trend can be con- 
sidered an argument against the dust disk clumping hypothesis. 
However, we cannot totally exclude that a random distribution 
of clumps may reproduce the observed visibilities. 





Parameter 


Value 


Po 


Photospheric density 


4xl0 -12 g.cnr J 


0(pole) 


Polar mass flux 


10-" M .yr' 


0(eq) 


Equatorial mass flux 


10" 10 M .yr -1 


mi 


Mass flux law exponent 


100 




Polar terminal velocity 


200 km.s -1 


Vco(eq) 


Equatorial terminal velocity 


0.1 km.s -1 


m 2 


Terminal velocity law exponent 


100 


i 


Inclination angle 


60° 


P.A. 


Major-axis orientation 


15° 


Teff 


Effective Temperature 


8250 K 


R* 


Stellar radius 


65 Rq 


d 


Distance 


650 pc 



Table 5. Parameter values for the SIMECA model. 



6.1.2. Inner gaseous disk 

To test if the inner gaseous disk of HD 62623 can be responsible 
for the residual visibility, we used the SIMECA code developed 
by Stee (1996) to model the circumstellar gaseous environment 
of active hot stars. This code is based on a latitudinal-dependent 
CAK wind model (Castor, Abbott, & Klein 1975) and large ve- 
locity gradient radiative transfer (Sobolev 1960) and is able to 
compute SED, hydrogen emission line profiles, and intensity 
maps in the continuum and in the emission lines. However, in our 
case we only need to compute radiative transfer in the continuum 
in order to obtain 8 and 13 pm intensity maps and 0.1-100 pm 
SED. The aim of this modeling is not to constrain the physical 
properties of HD 62623 's gaseous inner envelope, which is not 
possible due to our small number of measurements, but rather 
to determine if this component can be responsible for the 70 m 
N-band residual visibility. 

In the SIMECA code, mi and m 2 are exponent used to set 
the latitudinal dependence of the mass-flux (O) and terminal 
velocity (Vco), respectively: 

<D(6>) = ®(pole) + (®(eq) - <&{pole))sin m 8 (4) 
Voo(9) = Vco(pole) + (v^ieq) - v xl (pole))sin m2 8 (5) 

Typical values of the terminal velocity of an A super- 
giant wind are of the order of 500 kms -1 (Achmad et al. 
1997) although Lamer et al. (1995) derived a terminal ve- 
locity of 150 +50 kms -1 for HD 62623 from IUE observa- 
tions. However, HD 62623 may be rotating at a significant 
fraction of its breakup velocity; i.e., V sin i ~ 70kms -1 and 
160 kms -1 < v c < 250 kms -1 , depending on the value of M*. 
Thus, the expansion velocity at the equator can be smaller, as 
predicted by the bi-stability model (Lamers & Pauldrach 1991). 
Studying the classical Be stars a Arae and k CMa with the 
VLTI/ AMBER instrument, Meilland et al. (2007a, 2007b) found 
that equatorial expansion velocities were smaller than 10 kms -1 . 
Moreover, a spectroscopic follow-up of the variable Be stars 5 
Sco (Miroshnichenko et al. 2003) and Achernar (Kanaan et al. 
2007) leads to values of the order of 0.2-0.4 kms -1 . 

Thus, we have tested various geometries and kinematics to 
simulate HD62623's gaseous environment : a spherical wind 
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with terminal velocities between 10 and 500kms~', and a flat- 
tened envelope with terminal velocity between 100-500 kms" 1 
at the pole and 1-100 kms -1 at the equator. We have also tested 
various mass-loss rates and densities. This modeling provides a 
typical extension for the gaseous environment of 5-10 mas in 
the N-band (assuming a 65R star located at 650 pc). Thus, the 
visibility of such a component for a 70 m baseline is about 0.8- 
0.95 at 8 /im and 0.91-0.97 at 13 /.im. Finally, this inner gaseous 
disk should be almost non-resolved and can explain the observed 
residual visibility. 

Assuming that all other components are fully resolved, we 
can roughly deduce the 8 and 13 fj.m flux of this structure to be 
20% and 12% of the total flux, respectively. To produce such a 
contribution to the observed flux, the mass of the gaseous enve- 
lope has to be of the order of 4. 10~ 7 M . 

To illustrate this scenario we combine a MC3D and a 
SIMECA model without taking into account the interaction be- 
tween the gas and the dust. We calculate the complex visibility 
for each contribution separately and finally add both components 
using the formula : 



follows a power law as a function of the distance to the central 
star : 



ViFi + V 2 F 2 



(6) 



The resulting visibilities are plotted as a function of the spa- 
tial frequency in Fig [10] An 8 /mi SIMECA intensity map for 
the corresponding model is also plotted in this figure. The cho- 
sen density distribution is highly flattened in order to reproduce 
the visibility difference between B7 and Bg. To obtain a good fit 
of the interferometric measurements we need to slightly change 
the values of our MC3D best model since the compact gaseous 
emission increases the global visibility level. Thus, by increas- 
ing both the total mass to 3.5. 10 -7 M G and the inner dust ra- 
dius to 4.5 we managed to obtain a satisfying fit with a total 
reduce^ 2 for the nine baselines of 4.6. Table |5]present values of 
the SIMECA parameters corresponding to this model. Details on 
the physical meaning of these parameters are presented in Stee 
et al. (1996). Finally, Fig QT| presents a composite image of the 
MC3D + SIMECA model close to the central star. We note that 
the gaseous envelope is not strongly constrained and that these 
parameters have to be considered only as an example of a rela- 
tively "good" model. 

6.2. Physical conditions in the inner gaseous envelope 

One of the main issue regarding supergiants showing the B[e] 
phenomenon is how dust can form in the surroundings of these 
highly luminous, hot stars. To form dust in a circumstellar envi- 
ronment, two conditions must be met: the temperature must be 
lower than the dust sublimation temperature (i.e., between 800 
and 2000K, depending on the chemical composition of the dust), 
and the density of the chemical elements involved in the dust for- 
mation have to be higher than a critical value. 

These two conditions are hard to find at the same location 
in the circumstellar environment of these stars since, for most 
models, the density quickly falls with the distance. However, the 
stellar radiation might be shielded by the gaseous inner envelope. 
This would allow the temperature to decrease faster than in a 
normal reprocessing disk, and thus, the sublimation temperature 
could be reached closer to the central star. 

Since we managed to significantly constrain the inner rim 
extension of a supergiant dust disk for the first time, we can test 
this hypothesis by inferring the temperature variation in the inner 
gaseous envelope. We assume that the temperature distribution 



T(r) = T, 



eft'. 



R* 



(7) 



where T e ff an d R* are the stellar effective temperature and 
radius. In the vacuum, the only effect on the stellar radiation is 
the geometrical dilution of the energy and thus y-0.5. For a sim- 
ple reprocessing disc, i.e. stellar radiations absorbed by the disk 
and re-emitted at the LTE temperature, y-3/4- (Porter 2003). If 
the stellar radiation is highly shielded by the gaseous environ- 
ment, then y »3/4. 

Using R,„=3.85±0.7 AU (or 12.8+2.5*) and T,-„=1250K, we 
can determine the value of y for HD 62623 's circumstellar disk. 
Using equation (4) with r=R,„ leads to : 



fc) 



= 0.74 ±0.1 



(8) 



This value of y is very close to the value of the reprocess- 
ing disk. If we use the estimation of the inner radius from our 
SIMECA + MC3D model, i.e. R in =4.5 AU, we find y -0.79, 
again close to the reprocessing disk. Consequently, no strong 
shielding effect of the stellar radiation by the gaseous envelope 
is needed to obtain a temperature low enough to form dust at R,„. 
If the star is rotating close to its critical velocity, the equator can 
even be cooler thanks to the Von Zeipel effect. This would even 
decrease the value of the y parameter. 

6.3. Comparison with other A-type supergiants 

HD 62623 is one of the coolest stars that exhibits the B[e] phe- 
nomenon. Actually, it's a member of the very rare A[e] super- 
giant class. A similar object has recently been discovered in the 
Small Magelanic Cloud by Kraus et al. (2008). Apart from the 
spectral features relative to the B[e] phenomenon exhibiting the 
presence of a highly illuminated gas and dust region, these stars 
seem to share the same properties as their "normal" counterparts. 

Normal-A-Type supergiants are known to exhibit radiation- 
driven stellar wind with a terminal velocity of the order of sev- 
eral hundreds of kms -1 and mass-loss between 1O~ 8 M and 
1O~ 6 M (Achmad et al. 1997). However, no dust has been found 
in the circumstellar environment of these standard A-Type super- 
giants. As already discussed in the previous section, this implies 
that the critical density of chemical elements that can be involved 
in the dust formation (i.e., mainly Si, O, and C) is reached too 
close to the star, where the temperature is still higher than the 
sublimation limit. This is probably due to the fact that the rela- 
tively high wind terminal velocity is reached close to the central 
star. Aufdenberg et al. (2002) studied the mass-loss and radiative 
stellar wind of the A-supergiant prototype Deneb. They approx- 
imated the wind velocity assuming: 



v(r) = v x (l-R*/rf 



(9) 



and found that v M = 225 kms -1 and /3 = 3. 

However, even if Deneb and HD 62623 have the same 
spectral class (i.e., A2Ia), their physical parameters can 
be different enough to affect the wind properties. For in- 
stance, Schiller & Przybilla (2008) derived a luminosity of 
1.96xl0 5 L o for Deneb, thus about 12 times our estimation 
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Fig. 11. A composite map of HD 62623's circumstellar environment in the N band from our MC3D + SIMECA model. 



for HD 62623. Remembering that the mass loss rate scales 
the luminosity, this should imply to a stronger wind for 
Deneb. But, in their estimation the authors used a distance 
of 802 pc, while the one derived from "Hipparcos the New 
Reduction" (van Leeuwen, 2007) is 1.86 smaller (i.e., 432 + 
61 pc). Finally, using this latest estimation of the distance, we 
derived a luminosity of 5.5 ± 1 x 10 4 L , closer to HD 62623 
luminosity. 

Thus, we can still use Deneb's wind parameters as a 
rough estimate for HD 62623. Consequently, at 12.8 R*, i.e. 
dust sublimation radius of HD 62623 determined by our MIDI 
observation and MC3D modeling, the expansion velocity should 
be of the order of 176 kms -1 . 

It is difficult to explain why HD 62623 is forming dust 
in its circumstellar environment while Deneb exhibits a purely 
gaseous stellar wind. Among the putative explanations for such 
a difference is a possible effect due to the stellar rotation. In the 
bi-stability model developed by Lamers & Pauldrach (1991), 
the fast rotation induces a strong opacity change in the wind 
between the pole and the equator resulting in the existence of 
two distinct regions in the circumstellar environment; i.e., a 
quickly expanding, low density polar wind and a dense equa- 
torial outflow with a smaller terminal velocity. However, they 
conclude that this mechanism is more efficient for stars with 
15000K < T e f f < 30000K. For lower temperature, the mass- 
loss rate is usually too low to reach the optically thick wind near 
the equator unless the stars are almost critical rotators. 

In the case of HD 62623, considering that V sin i~ 70kms~', 
i=60+10°, and 160 kms _1 < v c < 250 kms^ 1 , the rotational ve- 
locity v ra , ranged between 0.3 and 0.6 v f , which is not large 
enough to reach the bi-stability limit and produce a dense, opti- 
cally thick equatorial region. Moreover, the efficiency of the bi- 
stability model also increases with the luminosity as shown 
by Pelupessy et al. (2000). For instance, for a 20000 K star 
rotating at 60% of its critical velocity, the polar-to-equatorial 
density ratio is 0.95 for L* = 10 4 5 L while it decrease down 



to 0.21 for L* = 10 6 L . Finally, other phenomena have to 
be present to explain the differences between a standard A-type 
supergiant and HD 62623. 

6.4. The putative binary nature of HD 62623 

Studing the periodical variations of HD 62623's radial veloc- 
ity (137.7 or 161.1 days), Plets et al. (1995) found evidence 
of the binarity of this star. They concluded that the system is 
composed of a massive primary; i.e., Mi=31M or Mi=39M , 
and a small mass ratio of M 2 /Mi=0.03 or M 2 /Mi=0.15. The 
estimation of the projected semi-major axis of the system is 
1.6 AU < a. sin i < 2.4 AU. Considering that the companion or- 
bits in the equatorial plan of HD 62623's primary, and using our 
estimation for the inclination angle i =60 + 10°, we obtain a 
semi-major axis of a = 2.3 ± 0.6 AU. 

Finally, Plets et al. (1995) proposed a general scheme for 
HD 62623. The matter is probably ejected by the evolved 
massive primary stellar wind with typical terminal velocity of 
80 kms _I . However, the matter located in the equatorial plane 
cannot reach this value since it is decelerated by the companion, 
finally forming a circumbinary disk. Thanks to this deceleration, 
the density in the equatorial plane could be 10 times higher than 
the typical density for a single A2I star stellar wind. 

Binarity can be detected using interferometric techniques, 
and recent VLTI/AMBER observations of the classical Be 
star 6 Cen by Meilland et al. (2008) and VLTI/AMBER + 
VLTI/MIDI observations of the unclassified B[e] HD 87643 by 
Millour et al. (2009) resulted in the discovery of companions 
around these two stars. However, HD 62623's companion can- 
not be detected in our VLTI/MIDI data since the N band flux 
is fully dominated by the dusty disk emission and because the 
projected separation between the two components of the system 
is too small compared to the spatial frequencies probed by the 
instrument. 
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Nevertheless, considering the stellar parameters determined 
by Plets et al. (1995), the luminosity of the central star is at least 
14 times larger than the companion's, since for Mi=31M and 
M2/Mi=0.15, we obtain M2=5M Q , and the companion spectral 
class is between B8III and B6 V. Thus, it would result in visibility 
oscillations with a maximum amplitude of 2/14 ~ 0.14 for the 
binary system assuming no envelope contribution. In the H and 
K band, the envelope emission should represent 20% and 60% 
of the total flux, respectively. Thus, the oscillations should have 
a maximum amplitude of 0.11 in the H band and 0.056 in the 
K band and could be detected, at least in the H band, with the 
VLTI/ AMBER. 

Finally, it is interesting to note that the putative companion is 
probably located inside the inner radius of the dusty circumstel- 
lar disk, as seen by the VLTI/MIDI instrument. The disk geom- 
etry determined using the MC3D code and our interferometric 
observations is thus compatible with Plets et al.'s (1995) scheme 
for HD 62623. 

7. Conclusion 

Thanks to these first VLTI/MIDI observations of HD 62623, we 
were not only able to determine the extension and flattening of 
a supergiant star showing the B[e] phenomenon, as was already 
done by Domiciano de Souza et al. (2007) on CPD -57° 2874, 
but also to significantly constrain the inner rim of the dusty disk 
and the disk opening angle. However, we were not able to infer 
the possible disk flaring. 

We have shown that the inner gaseous envelope likely acts as 
a reprocessing disk. Since we detect residual visibilities for the 
70 m baseline, it is also probable that the inner gaseous envelope 
free-free N-band emission is not negligible and contributes to 
10-20% of the total flux. However, other hypothesis like clump- 
ing within the disk cannot be totally ruled out. 

How dust can form around this highly luminous 
(i.e., L* ~ 17000 L ) star still remains unclear. Moreover, 
HD 62623 is one of the few known A[e] supergiants, and it is 
hard to understand what the differences are between this object 
and other stars with similar spectral classes, like Deneb, which 
only exhibits a moderate radiatively driven stellar wind and no 
dust in its circumstellar environment. Lamers & Pauldrach's 
(1997) bi-stability model does not seem to fit accurately 
HD 62623 data since this star is likely to have a moderate 
rotational velocity, i.e. 0.3<V e <0.6. 

Moreover, comparing the dusty disk inner rim extension 
R,„ = 3.8 ± 0.7 AU and Plets et al.'s (1995) companion orbit 
estimation; i.e., 1.6 AU < a.sin i < 2.4 AU, we can conclude 
that a stellar wind deceleration by the companion's gravitational 
effects remains the most favorable hypothesis. 

However, to definitely rule out other possible explanations, 
we need to constrain both the geometry and kinematics of the 
inner gaseous envelope. As already shown for classical Be stars 
(Meilland et al. 2007a, 2007b) and Herbig stars (Malbet et al. 
2007 ; Kraus et al. 2008), spectrally resolved interferometry 
with the VLTI/ AMBER instrument in medium (R= 1500) or high 
(R= 12000) spectral resolution mode is the most suitable method 
for that purpose. New VLTI/MIDI observations with baselines 
longer than 70 m can also bring additional constraints on the 
gaseous environment geometry. Finally, as shown for Be stars 
by Meilland et al. (2007a ; 2009) and Kervella et al. (2009), 
comparing N band VLTI/MIDI and H-K-band VLTI/ AMBER 
measurements can add additional constraints on the physical pa- 
rameters, for instance the temperature and density of the gaseous 
environment. 
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